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The regions of independent quantum states, maximally classically correlated states, and purely 
quantum entangled (supercorrelated) states described in a recent formulation of quantum 
information theory by Cerf and Adami are explored here numerically in the parameter space of 
the well-known exactly soluable Jaynes-Cummings model for equilibrium and nonequilibrium 
time-dependent ensembles. 
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The well-known exact solution of the Jaynes - Cummings (JC) model [1,2] describing 
the interaction of a single mode radiation field with a two-level atom is the basis for a vast array 
of the current experiments on foundations of quantum mechanics involving entangled states [3], 
new ideas for quantum optics [4], and novel device structures such as micromaser [5], single- 
atom laser etc. [6]. It is also a generic model of the interaction between two dissimilar quantum 
systems (Bose field and two-level-system, (TLS)) [7, 8]. Thus this model contains all the subtle 
features of quantum entanglement. A new quantum information theory for entangled systems of 
Cerf and Adami [9] (CA) mathematically classified (1) independent quantum states, (2) 
maximally classically correlated states, and (3) a new classically forbidden and therefore purely 
quantum entangled states called "supercorrelated" states. This is because they found that the 
conditional entropies of quantum entangled system can be negative, a feature missed in an earlier 
theory [10]. CA have subsequently explored many further consequences of this discovery [11]. 

The purpose of this paper is to employ the formal CA theory to elucidate these features 
numerically in the parameter space of ensembles constructed with the JC model, thus providing 
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an explicit display of the features that CA discovered. We use the exact eigen-solutions of the 
JC model for constructing two ensembles: (a) thermal (equilibrium) ensemble at a temperature T, 
where the JC Hamiltonian may be considered as a generic model (Bose-TLS) system and (b) a 
time-dependent ensemble of quantum optical interest [10, 12]. Two simple models of the 
radiation often used in the literature, (i) blackbody-like and (ii) pure coherent source [12, 4] are 
considered in (b). (For a discussion of ensembles, see [13]). For each of these models, we deduce 
the individual partial density matrices of the radiation and the atomic fields from the system 
density matrix. This provides us with explicit expressions to examine numerically the negative 
conditional entropy and information contained therein, as a function of the model parameters, a 
feature not considered in [9]. For these ensembles, we find the three regions of states described 
by CA. Unlike CA, we exhibit them in the parameter space of the model and the two ensembles, 
the temperature and time-scale of the system. 

In the notations of [2], the Hamiltonian of the JC model concerns a two-level atom 
(TLS) interacting with a single given mode of quantized radiation (boson) field of a given 
frequency, co, is described in terms of the usual creation, a , and destruction, a, operators of 
the boson field; the TLS is represented by the z-component of the Pauli matrix operator with the 
energy separation of the two atomic levels 0) o , and their mutual interaction is expressed in the 

rotating wave approximation: 

h a+r = h0)a a + ^^~°z +^ K \a o_ +a<J + J. (1) 

Here K is the dipole-interaction strength between the radiation and the atom or in the generic 

model, the boson-TLS interaction. Exact solutions of this interacting system are [1,2] 
H A+R\<P( n ' s )} = "r{0.(n,s)]\(p(n,s)), (s = 1,2) 

H A+R \0,g) = -h[(O o /2]\0,g), 

| (p(n, 1)) = Cos n \n + 1, g) + SinO n \ n, e), 

\(p(n,2)) = -SinO n \n + \,g) + CosO n \n,e). 

The entangled states {| s))}, (s — 1,2; n — 0,1,2....) and \0,g) are orthonormal and 

complete. Here s=l, 2 labels the entangled states which in their bare condition are the ground (g) 
and the excited (e) states, and n, the states of the boson field. Also 
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£l{n,s) = G)(n + y2) + {3-2s)X n , tan0„ = K^(n + l)/[(Aa))/2 + X n ], 
(Ad)) = (ft) - C0 o ), and A„ = ^(Aft)/2) 2 + K 2 (n + l). The angles 9 n are measures of the 
entanglement; 6 n — in the noninteracting case (no entanglement) and 6 n — 7T/4 for all n at 
resonance, Aft) = . For the sake of simplicity, the numerical results reported in this paper are 
for the resonant case. In Eq.(2), the eigenstates in terms of the complete set of states of the boson 
field and the TLS are also given. 

(a) Equilibrium ensemble: The density matrix constructed by using the maximum entropy 
principle subject to given average total energy of the system [2, 13], is thus found to be: 

oo 

Pa + r = E[|^(n,l))w(n,l)(^)(n,l)| + |^(n,2))w(n,2)(^(n,2)|] 

n=0 (3) 

+ \0,g)w(0){0,g\. 

Here w(0) = [exp(/?hft) /2)]/Z A+/? , w(n,s) = [exp- /Ml(n,s)]/Z A+R , <mdZ A+R is the 
system partition function. Here f3 — (k B T) 1 where T is the temperature. In terms of the 
boson and TLS states using Eq.(2), it is 

oo 

Pa + r = I \(\n + l>(n + l\)(\g)(g\)(w(n,l)Cos 2 e n + w(n,2)Sin 2 e n )] 

n=Q 

+ w(0)(|0}<0|)(|g)(g|) 

oo 

+ l\(\n){n\)(\e){e\)(w(n,l)Sin 2 6 n + w{n,2)Cos 2 6 n )\ 

n=0 

(4) 

OO 

+ ^(w(«,l) - w(n,2))Cos6 n Sin6 n 

n=0 

The third sum in the above represents the entanglement as well as the "decoherence" features of 
the interacting system. We obtain the "marginal" density matrix 

PA = Tr RpA+R = f(-)\g)(g\ + /( + )k)(4 (5) 

oo 

Here /(_)=w(0)+ ^yw(m,l)Cos 2 6 m +w(m,2)Sin 2 6 m j, andf^-l-f^ are the 

m=0 

occupation probabilities of the g and e states. Similarly, the "marginal" density matrix of the 
boson is 
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Pr = Tr ApA+R = (6) 

n=0 

Here p = w(0) + w(0,l)Sin 2 G + w(0,2)Cos 2 G , and for n=l,2,3, , 

2 2 2 * 2 

p n = w(n,l)Sin 6 n +w(n,2)Cos 9 n + w(n - 1,1) Cos n _x + w{n — \,2)Sin 9 n _\. 
From Eqs. (2, 5, 6), we compute S A+R = ~Tr(p A+R \np A+R ), S A = -Tr A (p A lnp A ), and 
5"^ = — Tr R [p R \np R ). Henceforth we focus on the resonant case CO — COq and the 
dimensionless variables of this ensemble are {/5hCo) 1 and (k/Co). 

(b) Non-equilibrium ensemble: We consider a unitary time-evolution (Liouville - von 
Neumann) of an initially prescribed density matrix of the system, p A+R (0), [1, 10, 12] given by 
PA +R (t) = U A+R (t)p A+R (0)U\ +R (t) where U A+R (t) = Qxp(-itH A+R /h). In terms of 

the eigenstates given in Eq. (2), 

r A+R 

(t) = \0,g)(0,g\p A+R (0)\0,g)(0,g\ + 

1 0, g) L (0, g \p A+R (0)| (p(m, s'))[exp it(co /2 + fl(m, s'))](p(m, s') | + /i.e. + 

£ | (p(n, s))((p(n, s) \p A+R (0)\ (p(m, /))[exp- if(£2(n, s) - ft(m, s'))](p(m, s') |. 

(7) 

For simplicity of presentation, we now employ two special initial density matrices, first 

considered in [12] and later used in [10] to examine interesting aspects of the entropy of the 

radiation field. There are two radiation models exhibiting different interesting behaviors [12, 4]. 

The initial state is specified in the form where the photons come from a single mode of radiation 

considered in [12] and the atom is in its excited (e) state: 

p A+R (0) = p A (0)®p R (0), p A (0) = \e)(e\, p R (0) = £ p(n)\ n)(n\ (8) 

n 

Model (i): Blackbody-like source 
p(n) - (N) n (l + N) ^ +n \ N - ^np(n) - mean number of photons. 

n 

Model (ii): A single mode pure coherent source [12] for which the photon number distribution is 

, , (NY exp(-Jv) / ^ 
Poissonian, p c {n) — ' v /n\' w * m ' anc ^ Pa\V) are as m model (i) . Note that 
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p c (n) has a peak at n — N whereas p(n) of model (i) is monotonic, but both vanish for large 



n. 



Using Eqs. (7, 8) we obtain 

[| <p(n ,\))p(n)Sin 2 d n (<p(n, 1) | + 1 (p(n, 2))p{n)Cos 2 6 n (<p(n, 2) | + 

PA+R(t) E 2|<p(nJ)){p(«)5m0„C^0 n C^[?(^(«a)-^(«,2))]}((p(«,2)^ ( ^ 

It should be emphasized that the preparation of the initial state as in Eq.(8) is nontrivial both 
conceptually and experimentally. However the above two models have served as guides and we 
use them here in the same spirit. In contrast to the equilibrium density matrix Eq. (3), the above 
expression contains only entanglement effects (dependent on K only). From this the "marginal" 



density matrices are 

Pa(0 = \g)w 8 (t)(g\ + \e)w e (t){e\, p R (t) = J»/> n (0(4 



(10) 



where w g (t) = Y,P( n W n (t), w e (t) = 1 - w g (t) and 

n 

P n (t) = P(n){\ - W n (t)) + p(n - l^^with W n (t) = (Sin 2 2d n )Sin 2 (tK^(n + l)). 
The total system entropy in both of these cases is time-independent and solely determined by the 



initial state density matrix, 

s A+R (t) = s A+R (t = o) = - 
s A+R (t) = s A+R (t = o) = - 



AnogN-(N + l)log(Af + l)] (model (i)) (11a) 



AHogAf - N -l> c («)ln(«!) 



(model (ii)) (lib) 



By construction this is just the entropy of the given initial radiation field. The counterparts of 

Eqs. (5, 6) are calculated using Eq. (10) for the two models. For the resonant case the 
dimensionless parameters for these two ensembles are| Kt 7Ta/TV j and N . 

The quantum conditional entropies express the residual information in the atomic and 

radiation systems respectively while retaining the quantum phase information are 

S(A + R\R) = S A+R -S R , S(A + R\A) = S A+R -S A . (12) 

As shown formally by CA, for quantum entangled subsystems these can be nonmonotonic and 
may even be negative unlike their counterparts in classical conditional entropies which are non- 
negative. And finally the quantum mutual entropy is S(A:R) — S A + S R — S A+R . As in the 
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classical case, S(A: R)>0 but can be <2mm[S A ,S R ]. When A and R are classically 
maximally correlated, the classical upper bound S(A:R) = min[5 , ^,5'^] is saturated. The 

range between the classical and quantum upper bounds corresponds to pure quantum 

entanglement and is called the state of supercorrelation [9]. This is the new feature found by CA 

and missed entirely in the earlier works [10]. In the present work, we will numerically explore 

these features as a function of the parameters of the model and the ensembles chosen, by 
focusing our attention on S(A + R\R)/S A for simplicity of presentation. Similar results obtain 

for S(A + R\A)/S R . 

Fig.l shows (S A+R — S R )/S A vs. (/5h&)) 1 for case (a) of the generic boson-TLS 
system for three representative values of (k/(o). For (fih(d) 1 — > the system approaches its 
ground state, thus (S A+R — S R )/S A — > 1 for all (k/G)), as expected. For (k/G))=0.5, 
(S A+R — S R )/S A remains positive whereas for (k/G))=2.5 and 5, negative regions (pure 

quantum entangled supercorrelated states) appear at finite temperatures. The region between the 
values 1 and zero of (S A+R — S R )/S A , represents maximally correlated states. 

For (k/Co) greater than 1, the interaction part of the Hamiltonian in Eq. (1) dominates 

over the TLS part. The eigenvalues £l(n, 2) then become negative for increasing values of n, as 

(k/Co) increases (for the examples considered here, Q.(n,2) is negative for (k/(o) =0.5, n=0; 

for (k/Q))=2.5, n=0, 1, ...5; and for (k/0))=5, n=0, 1, ...25) and so at low temperatures, the 

corresponding weights dominate the density matrix, reminiscent of chaotic behavior of the 
system as (k/Q)) increases. This is analogous to the result of Furuya et. al. [14], displaying 

chaotic features in a semiclassical version of the JC model, due to entanglement effects. If this 

result has to be applicable to the radiation- atom system, such strong dipole-interaction strength 

can only be achieved if the cavity-size in the quantum optics experiment is of nanometer 

dimensions and the atoms are cooled to very low temperatures as in [15]! 

In Fig. 2, we display (<S'^ +K — S R )/S A vs. ^Kt 7TV^V j for case (b), model (i) for three 

values of ./V. We observe that for N — l (as in the single atom laser case, [6]) 
{S A+R — S R )/S A is negative indicating the appearance of supercorrelated states for many 
values of [fCt Tl^N j whereas for = 50, (as for high mean photon number limit, [10]) 
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these supercorrelation states do not appear except in a very small region for low yKt 71^1 N J 
values. Indeed, (^ + ^ — S R )IS A for \Kt 71^[n\ near zero has the form, 



(the numerator is positive for the two radiation models considered here). This behavior is not 
perceptible in the Figures displayed. For N less than or equal to 5 the supercorrelated states 
occur but appear to be spread out; they are more spread out for ./V greater than 5, with a 
decrease in the number of supercorrelated regions as well as exhibit compression of the 
oscillations. There is no completely disentangled state as in case (a) by virtue of its construction. 

In Fig. 3, we similarly display the results for model (ii) of case(b), for the same three 
values of N . Here again N = 1 displays large number of supercorrelated regions whereas they 
decrease for N = 5. For N = 50, we have a very small negative region of (S A+R — S R )/S A 

near zero time-scale as mentioned above, but there is a new feature - oscillations (revivals) for 
longer times. Similar revivals were found and discussed in [3, 4, 10]. The positive regions of 



WA+fl — Sr)/Sa correspond to the maximally classically correlated states while we have no 

approach to the region of disentangled states, as expected. This is a manifestation of the 
difference between the two radiation models mentioned earlier. 

In conclusion, we have shown here by exploring numerically in the parameter space of 
the exactly solube JC model that the quantum supercorrelated behavior where the marginal 
entropy exhibits negative values in both the equilibrium (thermal) and nonequilibrium (time- 
dependent) ensembles in the corresponding parameter space, clearly demonstrating the 
entanglement effects contained in the JC model. We believe that by this explicit example, we 
have elucidated the formal results of CA concerning the existence of supercorrelated states in the 
quantum entangled systems. As noted here, for the JC model their presence depends on the 
ensemble considered and the interaction strength between the two subsystems. In view of the 
recent advances in single-atom quantum optics, especially as we approach nanometric cavity 
(e.g., quantum wells) sizes, we hope these fully quantum correlated states will be explored 
experimentally in the near future. 
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We wish to dedicate this paper to the memory of Professor Edwin T. Jaynes. Thanks are 
due to Dr. R. W. Rendell for drawing our attention to Ref.[13]. AKR and KLJ thank the Office 
of Naval Research for partial support of their work. 
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"An Entropic Description of Quantum Entanglement in the Jaynes-Cummings Model" Figure 1 

A. K. Rajagopal, K. L. Jensen, F. W. Cummings 

Fig. 1: (S^fl — S R )jS A vs. (/5h(o) 1 for three values of the scaled dimensionless coupling 
parameter, K/(0 — 0. 5, 2. 5, and 5. The inset exhibits the regions of cross over from positive 
to negative values of (S A+R — S R )/S A . 
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"An Entropic Description of Quantum Entanglement in the Jaynes-Cummings Model" 
A. K. Rajagopal, K. L. Jensen, F. W. Cummings 



Figure 2 



Fig. 2: (S A+R — S R )/S A vs. Kt 7t'\[N for blackbody-like source for three values of 

N = 1,5, and 50 . Note that the vertical scales are different for the three different cases. The 

insets exhibit the cross from negative to positive values of (S A+R — S R )/S A for small times. 
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"An Entropic Description of Quantum Entanglement in the Jaynes-Cummings Model" 
A. K. Rajagopal, K. L. Jensen, F. W. Cummings 



Figure 3 

Fig. 3: (<S'^ + ^ — S R )/S A vs. Kt 7Ta/ -/V for pure coherent source for three values of 
N = \,5, and 50 . Note that the vertical scales are different for the three different cases. The 
insets exhibit the cross over from negative to positive values of (S A+R — S R )/S A for small 
times. Unlike in Fig. 2, we have oscillations for N = 50. 
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